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High sensitivity receiver systems with near ideal polarization sensitivity are highly desirable for development of 
millimeter and sub-millimeter radio astronomy. Multimoded bolometers provide a unique solution to achieve 
such sensitivity, for which hundreds of single-mode sensors would otherwise be required. The primary concern 
in employing such multimoded sensors for polar imetery is the control of the polarization systematics. In this 
paper, we examine the angular- and polarization- dependent absorption pattern of a thin resistive grid or 
membrane, which models an absorber used for a multimoded bolometer. The result shows that a freestanding 
thin resistive absorber with a surface resistivity of 77 / 2 , where 77 is the impedance of free space, attains a beam 
pattern with equal E- and H -plane responses, leading to zero cross polarization. For a resistive-grid absorber, 
the condition is met when a pair of grids is positioned orthogonal to each other and both have a resistivity of 
77 / 2 . When a reflective backshort termination is employed to improve absorption efficiency, the cross-polar level 
can be suppressed below — 30 dB if acceptance angle of the sensor is limited to < 60°. The small cross-polar 
systematics have even-parity patterns and do not contaminate the measurements of odd-parity polarization 
patterns, for which many of recent instruments for cosmic microwave background are designed. Underlying 
symmetry that suppresses these cross-polar systematics is discussed in detail. The estimates and formalism 
provided in this paper offer key tools in the design consideration of the instruments using the multimoded 
polar imeters. 


1. Introduction 

Present astronomical instrumentation applications in 
the millimeter and sub-millimeter desire photon back- 
grounded limited sensitivities. There are two possible 
basic approaches to further improve the sensitivity by in- 
creasing the number of detected spatial modes received 
by an imaging system. The first is to build an array 
consisting of numerous single-mode sensors with high 
optical efficiency. The second, a multi-mode sensor, de- 
tects many spatial modes on a single sensor (see, e.g., 
Ref. [1] ) with well-defined angular and polarization char- 
acteristics. In astronomical observations, for example, 
both single-mode sensors and multimoded sensors (see, 
e.g., Refs. [2, 3]) have found wide use for radiometry and 
photometry. 

Various techniques can be used to specify the modes 
coupled to a sensor and the resulting system architec- 
tures can be categorized by their modal filtering tech- 
niques (Fig. 1). In the context of polarimetry at millime- 
ter wavelengths, where a significant use is for measuring 
the cosmic microwave background (CMB), the primary 
focus of the recent developments has been directed at 
large arrays of single-mode detectors with feed-coupled 
waveguide polarization diplexers [4-7] or planar-antenna 
coupled structures [8-13] which can be photolithograph- 
ically produced in large numbers. Here, polarization 
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diplexing is achieved on the detector chip for the hori- 
zontal and vertical single mode detector channels (Fig. 1, 
S3). Prior to these developments, dual-mode waveguide- 
based orthomode transducer (OMT) structures followed 
by single-mode detectors [14-16] were used to form a 
polarimeter following the traditions of microwave de- 
sign. Alternatively, single-mode dual-polarization sen- 
sors were realized by combining intrinsically multimoded 
polarization-sensitive bolometers (PSBs) with external 
modal filtering structures to define the angular accep- 
tance (Fig. 1, S2) [17-20]. Multimoded sensors for imag- 
ing have also achieved polarization sensitivity through a 
wire-grid analyzer [21-24] (Fig. 1, Ml and MP1). The 
analyzer grid architecture has also been employed in con- 
junction with a feed coupled array (e.g., Ref. [25]) to 
cleanly provide polarization sensitivity (Fig. 1, SI). 

A multimoded bolometric sensor that is intrinsically 
polarization selective, as opposed to those using the 
wire-grid analyzer, belongs to another class of sensors 
and opens up a new phase space for polarization sen- 
sitive instruments in millimeter wavelength (Fig. 1, M2 
and MP2). Such a bolometer would employ a thin re- 
sistive grid as a polarization selective absorber. An M2 
implementation is the sensor developed for PIXIE satel- 
lite mission proposal aimed at CMB polarization and 
frequency spectrum measurements [26] . This sensor em- 
ploys a pair of orthogonally positioned resistive grids, 
each of which is separately read out, attaining simul- 
taneous sensitivity to two linear polarizations. Other 
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Fig. 1. Various modal filtering techniques for polarization sensitive detector systems. Throughout this paper, a single-mode 
dual-polarization sensor denotes a system that consists of a pair of single-mode detectors that are sensitive to two orthogonal 
polarizations. SI: a single-mode dual-polarization sensor comprising two intrinsically-multimoded polarization-insensitive 
sensors. Here, single-mode limit is set by the feed-exit coupling port (symbolically shown as the diameter “d” of the port 
appropriately set compared to the wavelength) and polarization separation is achieved via quasi-optical grid diplexer. The 
blue structures surrounding each sensor correspond to cold baffling. S2: a single-mode dual-polarization sensor comprising 
intrinsically-multimoded polarization-sensitive bolometers (PSBs), where the single-mode limit is set in the same way as SI 
and the polarization separation is achieved via the patterned absorber structure of the bolometer. S3: a single-mode dual- 
polarization sensor where the radiation input is split into two polarization components by a single-mode orthomode transducer 
(OMT) and subsequently absorbed by a pair of detectors. The feedhorn here merely provides a controlled coupling of the single- 
mode set by the detector structure to a single-mode in free-space; it can be replaced by an immersion lens or a phased-array 
antenna structure. Ml: a multimoded polarimetor comprising polarization-insensitive multimoded sensors. The number-of- 
modes limit is set by the baffling limiting the solid angle seen from the sensor and the polarization separation is achieved via 
quasi-optical grid diplexer (the same as SI). MP1: a polarimetor array consisting of polarization-insensitive segmented sensors 
filling the space and a quasi-optical grid diplexer. The spatial-mode limit and polarization separation are achieved in the same 
way as Ml. Each pixel can be either single- or multi-moded depending on the size of the pixels relative to the wavelength and 
the solid-angle limit set by the baffling. M2: a multimoded polarimetor comprising a pair of multimoded polarization-sensitive 
bolometers, where the number-of-modes limit is set in the same way as Ml and the polarization separation is achieved via 
the patterned absorber structure of the bolometer. MP2: a polarimetor array consisting of polarization-sensitive filled-array 
sensors. The spatial-mode limit and polarization separation are achieved in the same way as M2. 


applications exploiting the unique features of this sen- 
sor are also proposed [27]. An MP2 implementation, 
an intrinsically polarization-selective sensor based on a 
pair of filled arrays, was explored in Ref. [28] where each 
detector absorber ~ A 2 on a side was realized from a 
patterned thin metal film on a Si membrane and spaced 
by ~ microns (Fig. 1, MP2). These types of sensors 
(M2 and MP2), similarly to Ml and MP1, require cold 
baffling to control the detector power loading, but share 
an improved mapping speed advantage compared to the 


single- mode sensors [29] . 

The focus of this paper is to provide an estimate of the 
angular- and polarization-dependent response pattern of 
the thin resistive absorber, and to show that bolometers 
corresponding to the M2 type in Fig. 1 can achieve low 
levels of polarization beam systematics. Further, the 
residual non-zero systematics are shown to have even- 
parity patterns. A potential use for these sensors is to 
probe the signature of the inflation in the early universe 
through the odd-parity, or so-called R-mode, patterns 
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in the CMB polarization [30, 31], which is not contami- 
nated by the even- parity beam systematics [32, 33]. Im- 
plications presented in this paper are general and ap- 
plicable to a wide range of devices using thin absorbing 
grids or membranes, even though our work is motivated 
by the specific implementation of the device mentioned 
above [26, 27]. We choose to analytically study this 
system for numerical investigations the approach and 
method of Ref. [34] could be adopted. The results we 
present would serve as a key tool in the design considera- 
tion for this class of sensors and the instruments employ- 
ing them. Our analysis gives attention to the response 
pattern at large off-axis incident angles. The use of the 
large off-axis angles makes stray-light control easier and 
is often desirable in maximizing the number of modes. 
The number of modes (N mo< j es ) is related to the absorp- 
tion area (A), the solid angle (12) and the wavelength 
(A) by AVI = IVmodes A 2 , and larger VL allows N modes 
to be maximized while restraining the increase of the 
physical detector size A. A combination of a fast final- 
focus lens and a multimoded polarimetor (Fig. 1, M2) 
would provide such a large solid angle, for example. We 
note that so-called hlled-array sensors (Fig. 1, MP1 or 
MP2) are often deployed in a similar configuration using 
a final-focus lens. For this reason, the results presented 
in this paper offer potentially useful considerations for 
the hlled-array sensors with thin resistive-membrane ab- 
sorbers [35-38] even though they may be used as dual 
polarization photometers in the single to several mode 
limit. 

There are a few assumptions in our analysis regard- 
ing the configuration of the device of interest. First, we 
assume the absorbing grid or membrane is resistive, as 
opposed to reactive, and the current in the direction nor- 
mal to the absorber plane can be ignored. These approx- 
imations are valid if the absorber’s physical thickness is 
small compared to the penetration depth of the resistive 
coating ( electrically thin , hereafter). Secondly, when the 
absorber is a pair of orthogonal resistive grids, which are 
sensitive to orthogonal polarizations, we assume that the 
pitch of each grid and the distance between the two grids 
are small compared to the wavelength, and thus they are 
effectively resistive sheets on the same plane. It can be 
shown that there is little near-field modal coupling be- 
tween the crossed grids in this regime. In Appendix D, 
we discuss the conditions on the physical dimension of 
the detector such that the above stated approximations 
are valid. For example, the device proposed in Ref. [26] 
satisfies such conditions. Under these assumptions, we 
can treat the grid pair and the membrane equivalently, 
except the grid pair may have different resistivity in the 
two orthogonal directions. 

We start by setting up a formalism to evaluate the po- 
larization systematics in Sec. 2. We review the standard 
measure of the polarization systematics such as cross 
polarization in a single-mode system, and extend them 
to a multimoded sensor that comprises an electrically 
thin resistive absorber. Section 3 provides a rigorous 


foundation of this extension through an S'-matrix for- 
malism. Notably, this discussion provides ground that 
generalizes the results presented in Sec. 4 to an arbi- 
trary incident mode, whereas the derivation in Sec. 4 
is based on plane-wave incident modes. In Sec. 4, we 
discuss the result of electromagnetic calculations for the 
response pattern of an electrically thin resistive grid or 
membrane absorber sheet with an infinitely large area, 
with and without reflective backshort termination. We 
then briefly discuss the systematics due to diffraction of 
a finite-sized absorber in Sec. 5. 

2. Polarization Systematics 

In the context of polarimetry, our focus in evaluating 
the goodness of a response pattern of a sensor is for con- 
trol of polarization systematics. In particular, low cross 
polarization and small differential polarization response 
are important for high polarization efficiency and low 
spurious polarization [32, 33], respectively. In this sec- 
tion, we describe the formalism to evaluate the ability 
of the sensor to reject these artifacts in the angular re- 
sponse and retain the polarization purity of the incident 
radiation. Throughout this section, the word sensor rep- 
resents the entire polarimeter including the coupling to 
the plane wave propagating in free space, while the word 
detector denotes a part that converts an incoming elec- 
tromagnetic wave to an electric signal that can be read 
out (Fig. 2). 

It is convenient to discuss the polarization system- 
atics of a multimoded system in contrast to that of a 
single-mode system, which has already been discussed 
in detail in literature. We will adopt Ludwig’s third def- 
inition [39], in which the electric field directions of the 
two linear polarization bases are 

e v — eg cos <f> — e<£ sin <j) , and 

e h = eg sin 0 cos <j > , 

where eg and are unit vectors in spherical coordi- 
nates (see Appendix A). As a convention, we define the 
z axis as the on-axis direction of the optics and take 
the standard definition of a spherical coordinate system 
specified by (r, 0 , (f>) where 6 is the angle from the z axis. 
The vertical (horizontal) polarization direction e v ( ) 
asymptotes to the x ( y ) axis direction for 9 — > 0. 

It is customary to characterize the beam of a single- 
mode sensor through the radiation process. Although 
relevant for a receiver instrument is the reception pro- 
cess, time-reversal symmetry allows us to naively relate 
the characteristics in the radiation process to those in 
the reception process. For a single-mode sensor nomi- 
nally sensitive to vertical [horizontal] polarization, the 
radiation field pattern is described using co- and cross- 
polar far-held functions, G V qq{0,^>) and G v ^p{9,<j>), as 

p —ikr c r , •* (2) 

{Geo ( 0 , <t>) e«[h] + Gxp(0, <£) e h [ v] j , 
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where k is the wave number (fc = 27 t/A) and the super- 
script v or h (vertical or horizontal) denotes the nominal 
polarization of the detector. In practical applications, a 
polarimeter is often a dual-polarization sensor equipped 
with two detectors nominally measuring two orthogo- 
nal polarizations. We refer to such a dual polarization 
sensor as a single-mode system in the limit each polar- 
ization sensitive detector couples a single mode and the 
polarization isolation can be treated as a subdominant 
perturbation to the sensor response. 

The beam systematics are quantified using the far- 
held functions. We define the cross-polar level XP(9) 
as cross-polar radiation power: 


XP(0) = max 

<t> 


G X p(9, 4>) 
G v co ( 0 , 0 ) 


( 3 ) 


where we indicate the worst case in <j> by max^. Note 
that we do not assume normalization of the far-held 
functions to the on-axis beam and thus the denomina- 
tor in Eq. (3). The lack of normalization, as opposed to 
what is often standard, is our intentional choice since it 
is more appropriate when describing a multimoded sys- 
tem as we see later. The differential response is dehned 
as the difference of co-polar radiation power between the 
two detectors: 


DR{9) = 


GhoiW) 

2 

G h co (0A) 

2 

G v co ( 0 , 0 ) 


G h co ( 0 , 0 ) 



( 4 ) 


where we again take the worst case in <f>. Zero system- 
atics, i.e., XP = 0 and DR = 0, are attained when the 
far-held complex gain response satishes 


G£ o (M) = Gk o (M) = G(0), 

G v xp{9 7 4>) = G XP (9, 4>) = 0 ■ 


We briefly comment on some general properties of the 
far-held functions and the systematics measures. These 
properties apply to a scalar feed as well as a multimoded 
absorber sheet with isotropic resistivity, which we will 
show later. If the sensor of interest is symmetric under 
a 90° rotation about the z axis, the far-held functions 
satisfy the following relation: 

= Gq 0 (9, (/) + 7t/2) 

G\p{e^) = -G h XP {e^ + tt / 2 ), 


since the 90° rotational operation acts on the elements 
in Eq. (2) as (E Vl E h ) (-E h ,E v ) and (e v ,e h ) ->■ 
(—e.h, e„). This justihes the dehnition of Eq. (3), which 
may otherwise appear different for p = h. When the 
coupling to free space possesses continuous rotational 
symmetry about the z axis, the <f> to achieve max^ in 
Eqs. (3) and (4) corresponds to <j> = 7r/4 and (j) = 0, 
respectively [40]. DR(9) corresponds to the difference 
of the co-polar beams for the E- and H-planes when this 
condition is achieved. It is known that a necessary and 


sufficient condition for a scalar feed to attain zero sys- 
tematics [satisfy Eq. (5)] is to posses symmetric E- and 
U-plane responses, or G v co {9, 0 ) = G v co (9,tt /2) [39-42]. 
The orders of magnitudes oi XP and DR are typically 
related as DR ~ \[XP since | G v co — Gqq\ ~ \G X p\- 

The radiation held pattern is widely used to describe 
the beam pattern and its systematics for a single-mode 
system. In contrast, in this paper we mainly discuss re- 
ception processes in deriving the beam pattern of a mul- 
timoded sensor. This is because a multimoded beam 
pattern can be represented as a combination of multi- 
ple radiation held patterns and the relative excitation 
strengths among the radiation modes, and the latter is 
more naturally derived through reception processes. We 
assume a multimoded sensor with an electrically thin 
resistive absorber (Fig. 2) lying on the x-y plane. In 
the reception process, an incident plane wave induces 
absorber-surface current j x ( j y ) in the x ( y ) direction. 
We relate the induced current, j p , and the amplitude of 
an incident plane wave, a p (0, (f> ), with an incident angle 
( 9 , (f>) and a polarization p by co- and cross-polar cou- 
pling coefficients, G p co [9,(j)) and G XP (9,cj>): 


( Pxjx \ = f G v co (6,<j)) G XP (9, 4>) \ f a v (9,4>) \ 

\ PvJv ) V G h XP {9 , ft G h co (9, ft a h{9, ft ) ’ 

( 7 ) 

where p x (y) is the absorber surface resistivity in the x(y) 
direction. One would build a dual-polarization sensor 
by using a stacked pair of resistive grids as the absorber 
sheet, where each of the grids is aligned in the x or y 
direction, and by coupling the thermal signals from j x 
and j y to different power sensors [26] . 

Since the system of interest is multimoded, it is im- 
plicit in Eq. (7) that the surface current has multiple 
excitation modes for each of j x and j y . The modes may 
be expressed in terms of two-dimensional Fourier modes 
of the surface current, and the modal content of the 
induced current is dependent on (9,<fi,p), the incident 
angle and polarization of the plane wave. We will clar- 
ify the details of the modal content in the next section. 
The coefficients G P co ^ xp ^[9, <j>) represent the couplings 
between the current and the plane waves regardless of 
the modal content. When the system of interest is single- 
moded, the modal content does not depend on (9,<p,p) 
and G P co y xp ^{9, <f>) is equated to the far-held functions 
G P co[xp] (0> •A) through the symmetry between the ra- 
diation and reception processes. This is not the case 
for a multimoded system; e.g., the G p co ^ xp ^{9,<j>) do 
not describe a radiation held pattern as in Eq. (2). As 
discussed in the next section, however, we can regard 
G P co[xp\ (^> < A ) equivalently to the far-held functions for 
an inhnitely large, electrically thin resistive absorber 
sheet. For example, we show below that one can substi- 
tute them into Eqs. (3) and (4) in place of G^O[ xp] (^> 0) 
to evaluate the beam systematics. 
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Fig. 2 . Schematic figures showing our definition of the detec- 
tor and the sensor, and the mode couplings to them, for the 
examples of (a) a single-mode feedhorn-coupled bolometer 
and (b) a multimoded bolometer using a sheet absorber. For 
the former, a feedhorn serves as a coupling to the electromag- 
netic waves in free space, and a bolometer serves as a detector 
that converts the electromagnetic waves to an electric signal. 
Here, a/i, a/2, •• • and 6/1, 6/2, • • • are the amplitudes of in- 
coming and radiating plane waves, respectively, and ad 1 and 
bdi are the amplitude of the radiation from the detector and 
that absorbed by the detector, respectively. An S matrix 
characterizes the coupling among these modes provided by 
the feedhorn. For a multimoded bolometer in (b), on the 
other hand, the absorber both provides the coupling to the 
free-space-propagating waves and converts electromagnetic 
waves to thermal signals. The thermal signals then propa- 
gate to the power sensor (e.g., thermistor or TES) and are 
read out as an electric signal. Panel (c) schematically models 
(b) and clarifies the two roles the absorber plays and the fact 
that multiple modes couples to a power sensor, where S in 
the figure implies the coupling among the modes. Here, the 
detector corresponds to the combination of the dissipating 
part in the absorber and the power sensor. 


3. Beam Characterization Using the Coupling Coef- 
ficients G 

In this section, we derive the relation between the cou- 
pling coefficients G p co , xp A9, <j>) and the beam charac- 
teristics, in particular the beam systematics, for an in- 
finitely large, electrically thin resistive absorber sheet. 
We first introduce a formalism using a scattering ma- 
trix, and relate it to the far-held functions of single- 
mode and multimoded systems. This allows us to gen- 
eralize the expressions of the beam pattern, and thus 
Eqs. ( 3 ) and ( 4 ), for a multimoded system. We then 
show in the limit of infinitely large, electrically thin re- 
sistive absorbers an incident plane wave with a specific 
(9,(f>,p) couples to only a single surface-current mode. 
This simplifies the generalized expressions back to those 
for a single-mode system, except the far-field functions 


Gp O [ XP ](0, </>) are replaced by the coupling coefficients 

Gp O [ XP ]( 0 , </>). This formalism also shows that our re- 
sult is applicable to an arbitrary incident mode, even 
though the definition of Eq. ( 7 ) and the drivation in 
Sec. 4 assume plane-wave incident modes. This gener- 
alization is of importance since spherical waves are the 
incident modes when a detector is placed at the focus of 
a telescope. However, the derivation in Sec. 4 does not 
require this groundwork. 


3. A. S'-Matrix Formalism 

A general response of a sensor can be described using 
a scattering matrix, or S matrix, relating the incom- 
ing and outgoing modes (see, e.g., Ref. [ 43 ] for a for- 
mal discussion) . Incoming modes are either plane waves 
from the sky propagating in free space, whose ampli- 
tudes are or radiation from the detector, 

whose amplitudes are a^i, a,d2, • • • . Outgoing modes, 
which are the time reversal of the incoming modes, are 
either plane waves reflected toward the sky with ampli- 
tudes 6/1, 6/2, • • • , or the modes absorbed by the detector 
with amplitudes bdi, &d2, ■ • • (Fig. 2 ). The amplitudes 
are proportional to the electric field strength and fol- 
low a Gaussian random distribution for radiation from a 
thermal source (e.g., the CMB). An S matrix describes 
the coupling of the incoming and outgoing modes pro- 
vided by an optical coupling element (e.g., a feedhorn) 
and relates the amplitudes as 

*>i{v) ='^2,S i:j {u) a j{u) , ( 8 ) 

i 

where the indices i and j runs over both /l, /2, • • • and 
dl,d2, For simplicity, we omit the dependence on 
frequency v hereafter. 

A single-mode sensor measures a single mode among 
the outgoing modes, which we label as dl. The S'-matrix 
elements Sdij relate the incoming plane- wave ampli- 
tudes a,j(j = /l,/2,---) and the measured amplitude 
bdi- 


bdi — 'y ^ Sdi j aj . 

The detected power P is calculated as 

P = (M 2 > = (I E S d ijaj\ 2 ) 

J=/b/2,- 

= E I^gI 2 <M 2 >> 

3=1 1 ,/ 2 ,- 


( 9 ) 


( 10 ) 


where the last equality holds when incoming plane wave 

amplitudes are uncorrelated to each other. We omit the 

coefficient, with the dimension of admittance, relating P 
■ 1 2 

and &di Since the index j runs over the plane wave 
modes with different incident angles and polarization, 

I 1 2 

Sdi j corresponds to the angular and polarization de- 
pendent antenna power pattern of the sensor in the re- 
ception process. 
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It is customary to define the far-field functions Geo 
and Gxp in the context of radiation patterns and thus 
they are related to the S-matrix elements for the radi- 
ation process, Sidi- The system of interest is usually 
symmetric under time-reversal operation [44] and thus 
the S matrix is Hermitian, Sij = ST. The Hermitian 
property of the S matrix allows us to equivalently re- 
late the far-held functions to the S-matrix elements for 
the reception process, S,j\ j , which are conceptually more 
straightforward to evaluate for a multimoded system. 

A natural extension of Eqs. (9) and (10) describes the 
response pattern of a multimoded sensor. The formalism 
differs from a single-mode system in that there are more 
than one mode coupled to the detector (Fig. 2). The 
detector-absorbed amplitudes bi are: 

bi= E Sij a j (* = dl,d2, ■ ■ ■ ) , (11) 

j=flJ2,- 


and cij — > a p (6,(ji), respectively. Equation (13) can be 
conveniently written in terms of a 2 x 2 matrix 


b= Jdfl S(9, <j))a(6, 0) , 


(14) 


with 


a{9,<t>) = 
S(M) 


a v (0,(/>) \ 

o-h (9, 4>) ) ’ 



( S™(M) s vh ( 0 ,<f>)\ 

\S hv (0,<l>) S hh (0,t) ) ’ 


(15) 


and dll = sin(9d(/>d(9. Note that the matrix S(9 , <f > ) does 
not have to be symmetric since it is an off-diagonal block 
of the entire S matrix, which is clear in Eq. (13). 

We then equate the matrix S(9,(f>) to the far-held 
functions: 


and the detected power P is: 


M E N 2 ) = ( E E S 


i=dl,d2,- 


ii CL n 


= E 

i—dl,d2,- 


E 


J=dl,d2,... _7=/l,/2,- 

2 

2 \ 


(M ) > 


(12) 


where, again, the last equality holds when the incoming 
plane wave amplitudes are uncorrelated. 

The absorbed- mode basis set, with amplitude bi ( i = 
dl, d2, ■ ■ ■ ), is chosen such that each of the modes cou- 
ples to an eigenmode of the dissipation process in the 
device (Fig. 2c). This corresponds to the order of op- 
erations in Eq. (12) that the square of the absolute 
value, |&i| 2 , is taken hrst and then summed over the 
index i. In the subsection 3.D we show that for an 
inhnitely-large, electrically-thin resistive absorber, two- 
dimensional Fourier modes of the surface currents serves 
as such eigenmodes. 


3.B. Far-Field Functions of Single-Mode Sensors 

Here, we relate the S matrix to the far-field functions 
of a single-mode dual-polarization sensor. For conve- 
nience, we label the two detectors and indices of the 
detector-coupled modes by their nominal polarizations, 
v (vertical) and h (horizontal), instead of dl and d2. 
Equation (9) is now rewritten as a pair of equations: 



Since the index j runs over free-space-propagating plane 
waves, we can choose the basis set such that the label 
j corresponds to a set of an incident angle (0, </>) and 
a linear polarization p (= v,h). Both S matrix and 
plane-wave amplitudes are relabeled as Sij -»• s* P (M) 


( S vv (6 ,</>) S vh (e,<t>) \_f G£ o (0,0) G™ P (M) \ 
\S hv (0,<t>) S hh (9, (j>) J \ Gx P (9, (j>) Gq 0 (9, </>) J ' 

(16) 

As noted earlier, G p co {9 1 <j>) and G P XP (9, (j)) slightly devi- 
ate from the standard definition in that they are not nor- 
malized to Gp O (0,0). Instead, their normalization car- 
ries additional information about absorption efficiency. 

Now we can rephrase the definitions of the beam sys- 
tematics in the context of the reception process. The 
cross-polar level, XP(9) defined by Eq. (3), is the re- 
sponse in power to a cross-polar incident wave normal- 
ized to the response to a co-polar on-axis incident wave. 
The differential response, DR(9) defined by Eq. (4), is 
the difference of the power response between the two 
detectors to their co-polar incident waves with the same 
incident angles. This language, in contrast to that of the 
radiation process, is directly applicable to a multimoded 
system, too. 

3.C. Far-Field Functions of Multimoded Sensors 

We relate the S matrix to the far-field functions and 
define the systematics measures for a multimoded dual- 
polarization sensor by extending those for a single-mode 
sensor. We relabel the index of the detector-coupled 
modes, dl,d2, ■■■, by a pair ( p,i '), where p (= v,h) 
labels each of the two detectors and the index i! runs over 
the modes coupled to each detector. Hereafter, we use 
symbols with primes for variables indexing the detector- 
coupled modes for clearity. 

As is done for the derivation of Eq. (14), we rewrite 
Eq. (11) as: 


U = f API {9 , <f))a(9 , 4>) , (17) 

J n 

where b‘ and S l (0, <j>) are defined similarly to Eq. (15) 
but with an additional index i' . We then define the 
far-field functions for a multimoded system as a natural 
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extension of Eq. (16): 

/ st(M) S* h (9,<i>) \ = f G v ^(e,<t>) G v j' p (0,<t>)\ 

\SL(W) ^(M) ) \G h jtUW) G h c l o(9A) ) 

(18) 

We note again that the normalizations of the far-held 
functions here are related to the absorption efficiency; 
they are not normalized to on-axis response. This is 
more important here than it is for a single-mode system 
since the relative normalization among G^O[yp] (^’ '?■’) 
for different i' tells us the relative excitation strengths 
among different detector-coupled modes. 

We adopt the definitions of the levels of cross polariza- 
tion and differential response rephrased in the context 
of reception process using response power: 


XP(9) = max ■ 


Ei' 




Ei' G&( 0,0) 


(19) 


and 


DR(0 ) = — max 
2 0 


E, G^(e,d>) 


Ei' G°cU 0,0) 


Ei- 

g££(M) 

CM 

Ei' 

Geo M 

2 


. (20) 

The only difference from a single-mode system is the 
use of Eqs. (12), (17) and (18), as opposed to Eqs. (10), 
(14) and (16), in calculating the response powers used 
in these definitions. 


3.D. Simplification for Infinitely Large, Electrically 
Thin Resistive Absorbers 

Here, we relate the far-held functions in a multimoded 
system G^, l o , xp J0,4>) with the coupling coefficients 

Gqo\xp](9,<P) i n Eq. (7) for an infinitely large thin re- 
sistive absorber, and show the equivalence between the 
G P CO [ X p\ (9, <t>) and the far-held functions for a single- 
mode system, G P co ^ xp j ( 9 , <f>). The surface current on the 
absorber sheet is denoted by its two-dimensional Fourier 
amplitudes, j x (k' x ,k y ) and j y (k' x ,k y ), where k! x and k' y 
are the wave numbers of the surface current, and j x and 
j y correspond to current density per unit area of ab- 
sorber in the x and y directions, respectively. For later 
convenience, we identify the Fourier modes using angu- 
lar variables (<?', <j>') through (k x ,k') = (kcoscj)', fcsin</>') 
and k = kcos9' (see Fig. 3); we relabel the current am- 
plitudes as j p {k' x , k' y ) — > j p (9', (j)') accordingly. We adopt 
the basis set of detector-coupled modes that maps di- 
rectly to the Fourier modes of the surface current and 
use {O' , (j)') as the index i' in Eq. (17): 




bvW,#) 

bh{9', (j>') 


Pxjx{9' ,(/)') \ 

Pyjy(9'A') J ’ 


(21) 




Fig. 3. Schematic figures showing the relation between 
the incident plane wave with a wave vector k = (k x ,k y ,k z ) 
and the induced surface current on an infinitely large ab- 
sorber sheet. For an on-axis incident plane wave (a), the 
electromagnetic field on the absorber surface is spatially 
uniform and thus the amplitude of the current is also uni- 
form, while the direction of the current is determined by 
the polarization of the incident wave. For an off-axis inci- 
dent wave (b), on the other hand, the surface-current am- 
plitude varies across the absorber sheet, which can be ex- 
pressed in terms of two-dimensional Fourier modes with wave 
vectors ( k' x ,k' y ) = ( k cos (j)' , k sin q 7) where (j)' is the polar- 
coordinate angle. Note that the amplitude propagation di- 
rection, (k x ,k y ), is independent of the direction of the cur- 
rent, ( jx,j y ); the former couples to the incident angle of the 
electromagnetic wave while the latter couples to its polariza- 
tion. For an infinitely large absorber, the wavelength of the 
current variation must be as long or longer than the wave- 
length of the incident wave (k < k), allowing us to express 
k as k = kcos6' . Further, the incident electromagnetic field 
on the x-y plane has a wave vector of (k x , k y ) and thus only 
couples to the current mode with ( k' x , k y ) = ( k x , k y ). In turn, 
this current mode only induces the electromagnetic field that 
has the wave vector ( k x , k y ). Thus, the coupling between the 
incident wave and the surface current is non- zero only for the 
Fourier mode of ( k x ,k y ) = ( k x ,k y ), or (8 ' = ( 8,(j > ). 


where p x and p y are the resistivity of the absorber in the 
x and y directions, respectively. We also relabel the S 
matrix, S l (9,<fi) S(6', <f>'\ 9, <j>). Recall that we omit 

possible frequency dependence and thus the discussion 
here is confined in a single wave number k. The de- 
tected power Pp for polarization p is calculated through 
Eq. (12) by substituting these definitions: 

Pp = J d Sl'pp \jp(9', cj)')\ 2 = -J dfi' \bp(9',<l > ')| 2 , 

(22) 

with dfl' = sin 9' dcj)' &9' and (p,p) = (v,x) or ( h,y ). We 
also rewrite Eq. (17) as 

b(9\ (j)') = I dH S(0', (j)'- 9 , </>)a(0, </>) . (23) 

Jn 

We comment on our choice of the basis set of the 
detector-coupled modes. In order for Eq. (12) to be 
valid, the basis set has to be the eigenmodes in dissi- 
pation process in the device of interest. A trivial choice 


here is the mode set where each mode couples to spa- 
tially localized surface current. Our choice to use the 
Fourier modes is equivalent to the trivial choice accord- 
ing to Parseval’s theorem (or the unitarity of Fourier 
transformation) . 

As shown in Fig. 3 and discussed in Sec. 4, an in- 
coming mode with an incident angle (0, <j>) only in- 
duces surface currents with a wave vector of (k' x ,k' y ) = 
(k cos 9 cos 4>,k cos 9 sin (f>). Thus, the S matrix can be 
written in a simplified form as 

S{9', 9 , $) = (5 (cos 9’ - cos 9) <S(0' - </>) S(9, </>) , (24) 


which is a product of a 2 x 2 matrix S that is only depen- 
dent on ( 9 , <j>), and the angular part of three-dimensional 
Dirac’s delta function in spherical coordinate system. 
Equation (23) now reduces to 


b(e,</>) = S(9, t j>)a(9,4>). (25) 

Substituting Eq. (21) into Eq. (25), we obtain 


f Pxjxifi) 0 ) 
V PyJyi M) 


S{9,4>) 


a v (9, (j>) \ 
o-h (9, <j>) ) 


(26) 


Comparing Eq. (26) with Eq. (7), one realizes the el- 
ements of the matrix S are equal to the coupling coeffi- 
cients previously defined: 


s(M) 


G£o(M) <? xp ( 0 , 4 >)\ 

G h XP (9A) Gq 0 (9, <j>) ) 


(27) 


The far-held functions G^q, xp -, ( 9 , (j>) are related to the 
S matrix via Eq. (18), where the index i' in the equa- 
tion is now labeled by (9',(f>'), and thus related to the 
coupling coefficients through Eqs. (24) and (27) as 


G CO[XP]( d,< ^ ( 28 ) 

= 5(cos9' - cos 9) 5(4> - cj>) G p CO y xp ^(9, </>) , 

As noted earlier, the G P co , xp , (0, </>) are not the far-held 
functions since they do not define a radiation held pat- 
tern by Eq. (2). The far- held functions are dehned by the 
left hand side of Eq. (28), where each mode labeled by 
(0', (f>') has a radiation pattern corresponding to a single 
plane wave, and G^ 0 ^ xp ^ (0, </>) represents the relative 
excitation strengths among the modes. 

Nevertheless, G P co ^ xp ^{9,(j>) hold a formal equiva- 
lence to the far-held functions in a single-mode system, 
Gcoixr] (^, ( t > )- For example, substituting Eqs. (25) and 
(27) into Eq. (22) leads to an expression for the detected 
power equivalent to that for a single-mode sensor: 


Pp — 


1 


dll 


Pp Jo 


G p co (9, 0)a p (0 , </>) + G P XP {9 , cf>)a p {9, 4>) , 

(29) 


with (p,p,p) = (x,v,h) or (y,h,v). Substituting 

Eq. (28) into Eqs. (19) and (20) yields forms of XP{9) 
and DR{9 ) that are equivalent to Eqs. (3) and (4) [45]. 
This formal equivalence to a single-mode system allows 
us to regard G P co , xp A9, </>) in Eq. (7), or in Eq. (28), 
equivalently to the far-held functions in a single-mode 
system. Note that Eq. (29) can be generalized to repre- 
sent the response to an arbritary incoming distribution, 
which can be expressed as a spectrum of plane waves. 

We note that our choice of identifying the Fourier 
modes by (9',<f>') has simplihed the derivation above, 
allowing us to evade possible Jacobians that would arise 
with an other choice, e.g., (k' x ,k' y ). Our choice is con- 
venient since, after relating (0',(j>') to (9,<j>) through 
Dirac’s delta function, the power measured by the de- 
tector can be calculated through an integral over the 
solid angle, P = JdSlp\j\ 2 , as is standard. This has 
allowed us to directly equate the matrix S(9 1 (f> ) with 
the G p ;o < xp i (9. (f>) defined in Eq. (7), which are inde- 
pendently introduced; the former is introduced in de- 
scribing the S matrix while the latter is the coupling 
coefficients in the process where a single plane wave in- 
teracts with the sheet absorber and induces the surface 
current. We also note that the derivation above has 
clarified the modal content of the current j p in Eq. (7) 
and its dependence on the incident angle (9,<j>)] in- 
duced current is a pure Fourier mode with a wave vector 
(. k' x ,k ' y ) = (fccos6»cos0,fccos0sin^). 

4. Response Pattern of an Infinitely Large, Electri- 
cally Thin Resistive Absorber 

In this section, we calculate the response patterns for an 
infinitely large, electrically thin resistive absorber. The 
absorber sheet may comprise either an electrically thin 
resistive membrane or a pair of orthogonal grids, where 
the latter could be constituted as a dual-polarization 
polarimeter by thermally coupling each of the grids to 
a power sensor. In the following, we treat the resistive 
grid-pair absorber equivalently to a membrane absorber 
regarding its electromagnetic coupling to free space, ex- 
cept that a grid-pair is more general in that it may have 
different resistivities in the two orthogonal directions. 
The equivalence is based on the assumptions that the 
grid is thin and finely pitched, and the distance between 
the stacked grids is small. Discussions on concrete con- 
ditions and the validity of the treatment can be found 
in Appendix D. Hereafter, we assume the more general 
case. We define the coordinates such that the absorber 
sheet is parallel to the x-y plane, where each absorber 
grid is aligned in the x or y direction, and the wave vec- 
tor of an incoming plane wave is k = (k x , k y , k z ) T with 
|fc| = k. Note that, although the calculation in this sec- 
tion assumes a plane-wave incident mode, the discussion 
in the previous section allows us to apply the calculated 
result to an arbitrary incident mode through Eq. (29). 
This is of importance since, in most of uses, the sensor 
is at the focus of optics, where the mode incident on the 
sensor is not a plane wave. Appendix E also discusses a 
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Absorber 




Fig. 4. The geometric arrangement of the incident, reflected, 
and transmitted fields and the thin resistive absober sheet 
placed at z = 0. 


case where the incident mode is not a plane wave. 

In the following, we first discuss the case of a free- 
standing thin absorber. Although this is a configuration 
that minimizes the polarization systematics, the max- 
imum absorption efficiency is limited to 50% as in an 
analogous case of the transmission line (see also Fig. 13 
in Appendix B). We then discuss the case that employs a 
reflective backshort, which is a technique frequently used 
to improve the absorption efficiency (see also Fig. 14 in 
Appendix C). As we will show, it significantly improves 
the efficiency while introducing small systematics. 



o 30 60 90 

9 (degrees) 

Fig. 5. Cross-polar level (top) and differential response (bot- 
tom) of an electrically thin, infinitely large resistive absorber 
with p = p, 0.7p, and 0.55?). Both are zero when p = p/2. 
One can see the approximate order of magnitude relation of 
DR ~ y/XP holds here. 


with 


4. A. Freestanding Electrically Thin Absorber 

We consider an arrangement shown in Fig. 4, where an 
electrically thin resistive absorber is at z = 0, lying on 
the x-y plane. The electric field of an incoming plane 
wave with arbitrary polarization can be expressed as 
a linear combination of vertically and holizontally po- 
larized components with amplitudes a v = a v (9,(j)) and 
= a/i ( 9 , (f i) . 

Ei(r; k) = e ~ lk r ( a v e v + a h e h ) . (30) 

The transmitted wave (E t ) and reflected wave (E r ) have 
wave vectors of k and k = (k x , k y , ±k z ) T . At z = 0, all 
three plane wave components have the same phase of 
e -i(k x x+k y y) anc j Qn iy ^.j ie sur f ace current mode with a 
wave vector of (k' x ,k' y ) = ( k x ,k y ) is induced. Thus, the 
current density of the induced current is written as 

J (x, y , z) = S{z)e~ l{kxX+k y y) (j x e x + j y e y ) , (31) 

where j x and j y are the Fourier amplitudes of the surface 
current per unit absorber area. In turn, this current only 
induces radiating electromagnetic waves with the phase 
0 f e -i(k x x+k v y) a t ^ = o (Appendix B), and thus the S 
functions in Eq. (24) are confirmed. 

Solving boundary conditions (see Appendix B), we de- 
rive the coupling coefficients in Eq. (7) relating ( j x ,j y ) 
and (a v ,a h ): 

( G v co /px G\ P /p x \ 1 = 1 / a+(p x ) P{p y ) \ 

\G h XP /p y G h co /p v ) 4 COS e\P(Px) a -(p y ) ) 

(32) 


a±(p ) = (2 p + rj)( 1 + cos 6) ± cos 2 ^)( 2 p — 77 ) (1 — cosf?) , 
f3(p) = sm2(j)(2p — 77) (1 — cos0) , 

(33) 


where p is the impedance of free space. When the 
resistivities are the same for the x and y directions, 
Px = P y = P, the cross-polar level and the differential 
response defined by Eqs. (3) and (4) are 


XP[9) 

= 7 cos 2 9(1 
4 


cos 9) 2 


(2p-p) 2 (2p + p) 2 

(2 p + pcos9) 2 (2pcos6 + rj) 2 

(34) 


and 


DR(6) = - sin 2 29 


\2p-p\(2p + p) 3 


(2 p + pcos9) 2 (2p cos 9 + rf) 2 


(35) 


Figure 5 shows XP(9) and DR(9) for various p. 

The most important implication is the case with p x = 
p y = p/2. As can be seen from Eqs. (32), (33), (26), and 
(27), this is the necessary and sufficient condition for the 
matrix relating ( j x ,j y ) and (a v ,a,h) to be diagonal, 


( jx A 2 V lcosS f a v \ 

\jyj 1 + COS 9 \ a h J 


(36) 


and for Eq. (5) to be satisfied for G/~, 0 ^ xp ^(9, (/>). Thus, 
this surface resistivity achieves zero cross polarization 
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and zero differential response, XP(6) = DR(9 ) = 0. We 
discuss the underlying symmetry that causes the sys- 
tematics to vanish later in this section. 

In this special case of p x = p y = 77/2, the angular 
dependence of the power absorption per unit area of the 
resistive surface, P a , is 

= = = (37) 


where ( p,j,a ) represents either ( p x ,jx,a v ) or ( p y ,jy,a h )• 
The power flow of the incident plane wave per unit area 
per polarization is Pi = r)~ 1 a 2 . Thus, the angular de- 
pendence of the absorption efficiency is 


P a {0) _ 2 cosfl 

cos 9 Pi (l + cos0) 2 ’ 


(38) 


where the extra factor of cos 9 arises in converting the 
power absorption P a (9) to be per unit area of incident 
wave. The efficiency is maximum at 50% at 9 = 0, as 
expected from an analogous transmission line configura- 
tion (see Fig. 13 in Appendix B), and slowly drops as 
9 increases since the effective surface resistivity for an 
off-axis incident wave deviates from 77/2. The antenna 
reception power pattern P{9) is obtained by normalizing 
P a (9) to its maximum: 

P( 0 ) = = 4cos2g 730) 

U "ma x 0 P a {9) (1+cos 0) 2 1 1 j 

where maxgE a (0) = P a ( 0) in this case. Figure 6 (top) 
shows the efficiency e(9) and the antenna power pattern 
P(9). Since e(9) is nearly constant, P(9) is close to cos 9, 
the naive expectation due to the geometric factor. By 
integrating P(9), we obtain the solid angle as a function 
of the maximum acceptance angle # max : 


fl(9 n 


0 =27T 
= 47T 


/ d(cos9)P(9) 

'0 

(3 — F 2 COS Umax) ( 1 COS 9 n 


1 + cos 9 max 

1 

+ 4 log - 


(40) 


- COSC 


Figure 6 (bottom) shows fl(0 ma x) in comparison to the 
naive geometric expectation of 7r sin 2 0 max . For a max- 
imum acceptance at the 9 max = 7 t/ 2, fl ~ 0.917T while 
the naive geometric expectation yields 7r. 

It may appear non-trivial that the response of an 
electrically thin resistive grid or membrane completely 
eliminates polarization systematics when p = 77/2. One 
can intuitively understand why the response pattern at- 
tains symmetry between the E- and E-planes (leading 
to zero systematics) as follows. As shown in Fig. 7, an 
E-plane (E-plane) incident wave is a linearly-polarized 
plane wave with its electric field parallel (perpendicular) 
to the plane of incidence. When p x = p y = p, the E- 
and E-plane coupling coefficients, G v co (9,(j) = 0) and 
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Fig. 6. Top: the absorption efficiency e(9) and the antenna 
power pattern P(9) of a large electrically thin resistive ab- 
sorber with surface resistivity p = 77/ 2 [see Eqs. (39) and 
(38)]. Also shown for comparison is the antenna power pat- 
tern for constant absorption efficiency, cos 6, naively expected 
from geometry. Bottom: the solid angle as a function 
of maximum acceptance 9 m ax for the thin resistive absorber 
with p = 77/2. Compared is the naive geometric expectation 
of 7r sin 2 0 max . 


Gp O (0,() = 0), relate the amplitudes of the incident 
field and the induced current as 

a v = — (p + 77 cos j x , 

cos 9 V 2 / / a 1 \ 

If n T]\ . 

ah = - a ( p cos 0 + — j v , 

cos0 V 2/ 

where a v (ah) is the electric field amplitude for .E-plane 
(E-plane) incident since e v = eg (eh = e^) at <j> = 0. 
One can see each of Eq. (41) consists of two terms. The 
first term dominates in the limit p 77 (thin film resistor 
approximates an open) and the second dominates in the 
limit p <C 77 (termination approximates a short). At the 
limit of p 77 (Fig. 8a), the electromotive force and 
current-induced field can be ignored. In this limit, the 
current is simply proportional to the component of the 
electric field parallel to the resistive surface, and written 
as 


p\J\ = cos9 E \E l \ (42) 

where 9e is the angle between the electric field and the 
absorber surface. Figure 7 defines 9e such that 9e = 
9 (9e = 0) for E-plane (E-plane) incident, and thus 
Eq. (42) is in agreement with the first terms of Eq. (41). 
At the limit of p <C 77 (Fig. 8c), on the other hand, the 
absorber approximates a perfect electric conductor, and 
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Fig. 7. Schematic figure of (a) .E-plane and (b) H- plane 
incident waves with wave vector k. The former (latter) has 
its electric (magnetic) field parallel to the plane of incidence. 
Shown are the relation among the incident angle (9), the 
angle between E field and the absorber plane ( 9e ), and the 
angle between H field and the absorber plane ( 9h )• 


the field of the incident wave and the current are related 
as J = 2n x Hi, where n = —e z is the unit vector 
normal to the absorber sheet. Thus, 

^\J\=cos0 H \E i \ (43) 

where rj\Hi\ = \Ei\ and 9h is the angle between the 
magnetic field and the absorber sheet surface. For FI- 
plane (FI-plane) incident, Oh = 0 ( Oh = 0) as shown in 
Fig. 7. When the value of p is at neither of the limits, 
the electric field simply equals to the sum of the voltage 
drop due to resistivity, Eq. (42), and the electromotive 
force, Eq. (43); the electric field of the incident wave is 
related to the current as 

(44) 

This follows from the boundary condition in forming 
the solution as a linear combination of the limited cases 
shown in Fig. 8a and 8c. The total incident field is the 
sum of those for Fig. 8a and 8c: Ei = E a + E' h ,. For 
this linear combination, the electric field projected on 
the absorber sheet is simply that of E a since it vanishes 
for Fig. 8c. Thus, Ohm’s law corresponds to E a \ s = pj , 
where | denotes a projection on the plane of the ab- 
sorber sheet. On the other hand, the field in Fig. 8c 
satisfies 2 n x H\, = J , where H\, is the magnetic field 
associated to E\, . Substituting these two conditions into 
Ei = E a + E l b , leads to Eq. (44) . 

Equation (44) is equivalent to Eq. (41). Now it is 
obvious why the current response to the incident wave 
becomes symmetric between E- and FZ-planes when 
p = 77/2, since this condition makes Eq. (44) symmet- 
ric to the exchange of E and H. This also clarifies the 
observation and can be used to understand why a wire- 
grid made from perfect conductor is not a perfect polar- 
izer [46]. 

It is worth noting that the same angular response as 
that obtained in the limit of p rj would also show 
up if we replaced the resistive sheet by a perfect mag- 



Fig. 8. Two ways to form a solution for a resistive absorber 
sheet interacting with an incident plane wave. One is to 
form a linear combination of (a) and (b). Here, (a) is the 
plane wave propagating without any interaction, or the limit 
of highly resistive absorber (p p), and (b) is the current 
on the sheet and the current-induced radiating fields. The 
boundary condition, or Ohm’s law, constrains the coefficients 
of the linear combination by ( E a + E b ) | = pJ, where E | 

denotes the electric field projected on the plane of the ab- 
sorber sheet. Another way to form the solution is to use a 
combination of (a) and (c), instead of the combination of (a) 
and (b). Here, (c) corresponds to a solution for the case when 
the sheet is a perfect conductor, as described in the text. A 
linear combination of (a) and (c) can form the solution since 
(c) is a linear combination of (a) and (b). 


netic conductor. In this sense, zero systematics are at- 
tained when there are both electric and magnetic con- 
ductor sheets and the contributions from the two are 
equal. This is exactly the same symmetry as the one 
discussed by Koffman [47] in the context of cross po- 
larization of a feedhorn, where the feedhorn coupled 
to a paraboloid reflector attains zero cross polarization 
when the radiation pattern can be described as a bal- 
anced sum of electric-dipole and magnetic-dipole radi- 
ation fields (Huygens source); see Appendix E for de- 
tails. The symmetry between the electric and magnetic 
dipoles discussed by Koffman corresponds to that of elec- 
tric and magnetic conductor sheets discussed here. We 
also note that our work can be seen as an extension of the 
current sheet model for phased-array antennas [48, 49]. 
The model assumes current sheet of surface impedance 
77 backed with a magnetic reflector to simulate a phased 
array, and derives the symmetry of antenna properties, 
such as the reflection coefficent and scanning impedance, 
between E- and IF-planes. It parallels to our discussion 
above and its result is consistent with our derivation. For 
example, it suggests that the absolute value of the reflec- 
tion coefficient is tan 2 | in both E- and IF-planes, while 
Eq. (38) can be rewritten as e(6) = (l — tan 2 |) /2. The 
factor two differences both in the absorption efficiency 
and in the resistivity to achieve the symmetry follow the 
expectation due to the absence of the magnetic reflector 
in our setup; see our discussion at the beginning of this 
section pointing out the analogy to the case of transmis- 
sion line (cf. Fig. 13 in Appendix B) and the maximum 
absorption efficiency of 50% without a backshort. 
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^ Absorber Backshort Absorber Image of absorber 




Fig. 9. (a) The geometric arrangement of the incident and 
reflected fields, the absober sheet placed at z = —d, and a 
reflective backshort placed at z = 0. (b) The virtual ar- 

rangement exploiting the image theory, equivalent to (a) in 
the region of 2 < 0. 


4.B. Electrically Thin Absorber with Reflective 
Backshort Termination 

Use of a reflective backshort behind the electrically thin 
absorber sheet improves the absorption efficiency and is 
desirable in practical applications. We consider an ar- 
rangement shown in Fig. 9a, where the absorber is at 
z = —d and a reflective backshort is at z = 0, with both 
parallel to the x-y plane. In calculating the response 
pattern, we exploit the image theory and use the ar- 
rangement of Fig. 9b. Solving boundary conditions (see 
Appendix C), the coupling coefficients are derived as 

{ Gco/Px Gjcp/Px \ 

\Gxp/Py Gqo/Pv ) 

= z~ ik ‘ d ( a+(p*) P hs (Py) \ 

4sinfc^dcos0 \ /3 hs (p x ) a°J{p y ) ) ’ 

(45) 

with 

a bs (p) = [(p + p) sin k z d — ipcosk z d] (1 + cos0) 

± cos 20 6 hs (p, d, 6) , (46) 

P hs (p) =sin2 <j> 6 hs (p , d,9) , 

and 

<5 bs (p, d, 9) = [(p — p) sin k z d — i p cos k z d\ (1 — cos 9) . 

(47) 

This equation is already instructive even without calcu- 
lating the cross-polar level and the differential response 
explicitly. In order to eliminate these systematics, the 
matrix elements should satisfy a5f = a bs and /3 bs = 0. 
This is only possible if <5 bs (p, d, 9) = 0. Even though the 
first term in 5 hs (p,d,9) vanishes for p = p, 5 hs (p,d,9) 
cannot be zero for an arbitrary 9 since k z = kcos9. 
Thus, there always are non-zero systematics, unlike the 
case without a backshort. However, by optimizing the 
distance to the backshort d , a low level of beam system- 
atics can still be achieved. 

For p x = p y = p, the cross-polar level and the differ- 
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Fig. 10. Cross-polar level (top) and differential response 
(bottom) of electrically thin resistive sheet absorption pat- 
tern with reflective backshort, for d = A/4, d = 1.2A/4, 
d = 1.4A/4, and d = 1.6A/4, where A and d are the wave- 
length and the distance between the sheet and the reflective 
backshort. 


ential response are 

XP(9) = \ K ^ n cos 2 p (1 _ cog p) 2 
4 sin kd 

Ipsinfcd — ipe lkd \ 2 \rjsmk z d + ipe lkzd | 2 

|p sin k z d cos 9 — ipe lkzd \ |psin k z d — ipe lkzd cos 9 1 2 

(48) 


o sm ka 

| p 2 sin 2 k,d — p 2 | |psinfcd — ipe lkd | 2 

x : 1 : — ; 2 ; - 1 2- 

|p sin k z d cos 9 — ipe lkzd \ |p sin k z d — ipe lkzd cos 9 \ 

(49) 

Note that Eqs. (48) and (49) are equivalent to the no- 
backshort case [Eqs. (34) and (35)] only if one sub- 
stitutes k z d = kd = tt/2 and p — ► 2 p. Figure 10 
shows XP(9) and DR{9) for p = p and various d. For 
d ~ 1.4A/4, a low level of the beam systematics, XP 
below — 30 dB and DR below — 20 dB, can be achieved 
for a wide range of acceptance angle, 9 < 55°. The XP 
and DR do not significantly degrade for d = 1.2A/4 or 
1.6A/4, suggesting a fractional bandwidth of ~ 30% is 
possible while maintaining the low levels of systematics. 

We define the angular dependence of the power ab- 
sorption per unit area of the resistive surface, P a , as an 
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average over <fi: 


P a (d) = ±-J^dcj>P a (9,<j>). (50) 

For one grid running in x direction, P a (9, (j>) here is 


P a {0,4>) = PxJx 2 


1 

V 
1 

V 


G v co (9, (f>)a v + G XP (9, 4>)dh 


Gco(9,(/)) + G XP (9,(f>) 


(51) 


where the last equality assumes unpolarized incident 
light with (|a„| 2 ) = (|a/i| 2 ) = a 2 and (a*dh) = 0. When 
Px = Py = V: Pa{0) IS 


P a (9) = sin 2 k z d cos 2 9 p 1 a 2 

x [4 sin 2 k z d( 1 + cos 9) 2 + cos 2 k z d (3 + cos 2 9)] 

1 

X t 

sin k z d (1 + cos 9) 2 + cos 2 k z d 

1 

X t ■ 

sin k z d (1 + cos 9) 2 + cos 2 k z d cos 2 9 

(52) 


Again, Eq. (52) is equivalent to Eq. (37) if one formally 
substitutes k z d = 7t/2 , except for an overall factor of 
two. Substituting Eq. (52) into the first equalities of 
Eqs. (38) and (39), we define the absorption efficiency 
e(9) and the antenna power pattern P(9), respectively. 
Similarly, we define the solid angle fl(0 max ) by the first 
equality of Eq. (40). Figure 11 shows them for vari- 
ous backshort position d for p x = p y = rj. Unlike the 
case without a backshort, 9 = 0 does not always give 
the maximum P a (9); e.g., P a ( 0) 0 ma xq P a (9) when 
1.4A/4 < d < A/2. 

In summary, the use of the reflective backshort leads 
to significant improvement of the absorption efficiency 
and maintains a large sensor acceptance angle. On the 
other hand, it also leads to non-zero cross polarization 
and differential response. However, these systematics 
can be suppressed for a wide range of incident angles, 
0 < 9 < 60°, by matching resistivity to the free-space 
impedance and optimizing the backshort position. We 
also point out that Eq. (46) shows that the differential 
response (sourced by the d bs term in a bs ) depends on </> 
as cos 2(j>, yielding an spurious polarization pattern with 
even parity only and maintaining a systematic-free mea- 
surement of the odd-parity patterns [32, 33]. This is of 
importance in the context of CMB polarization observa- 
tion, where the precision measurement of the so-called 
S-mode, or parity odd, pattern is the primary goal. 

In order to see why the beam systematics are small 
and yet what prevents their exact elimination, we look 
into the beam patterns in the E and H planes. For 
0 = 0, G h XP (9 , 0) = G V XP (9 , 0) = 0. When p x = p v = p, 
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Fig. 11. Absorption efficiency (top), antenna power pattern 
(middle), and the solid angle (bottom) for with p x = p y = r/ 
and with a reflective backshort. Cases with the backshort 
positions of d = A/4, d = 1.2A/4, d = 1.4A/4, and d = 1.6A/4 
are shown. 


the coupling coefficients are 


Gco(9i 0 ) 


G h co (9, 0) 


e ikzd 2p cos 9 

p cos 9 — ipe~ ik * d (sin fc 2 d) _1 
e lkzd 2p cos 9 

r] — ipe lkzd cos 9 (sin k z d) 1 


(53) 


We reform them and relate the incident field and the 
induced current using Eq. (7) as: 


CLh = 


g—ik z d 

2 cos 9 

g—ik z d 

2 cos 9 


p cos 9 + 
V 


1 — e~ lkzd cos k z d 
pcos9 


1 — e~ lkzd cos k~d 


(54) 


3v 


where a v (ah) is the electric- field amplitude for E-plane 
(H-plane) incident since <j> = 0. 

Equation (54) has an almost equivalent form to the no- 
backshort case of Eq. (41) except for an extra factor of 
—e~ lkzd cos k z d in the second term; without this factor, 
the beam pattern is symmetric between E and H planes 
when p = rj. As we see below, this extra factor can 
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with p = r] would recover the symmetry between E and 
H planes, or the symmetry between Eqs (42) and (43), 
if there were contributions only from (a) and (c). The 
relation of the held and current in (c) is the same as that 
of Eq. (42), except the current here is twice as large since 
E r a , is phase coherent to E l a at z = —d: 

p\J a >\=coB0E(\E i a \ + \E r a ,\)=2cos9 E \E i \ . (55) 


(c) (d) 




Fig. 12. Schematic figures of the field interacting with an 
electrically thin absorber and a backshort, (a) at the limit of 
p -C rj and (b) at the limit of p p. The latter is further 
decomposed into two: (c) where the total field is the sum of 
the incident field E l a and a component of the reflected field 
E' a , that is phase coherent to the incident field on the surface 
of 2 = —d, and (d) where the field is a superposition of E\,, 
(= -K') and El„ (= E r a ). 


be understood as a reception pattern of a two-element 
phased-array antenna. Note that the overall phase of 
e -ik z d gq. (54) is physically irrelevant, being due to 
our convention here that the resistive sheet is placed at 
2 = —d, not at z = 0. 

We again look into the terms in Eq. (54) by consid- 
ering the limits of p <C p (approximately shorted) and 
p p (approximately open). Figure 12 illustrates these 
limits. The first term, which dominates in the limit of 
p -C rj, corresponds to a reflection off a perfect conduc- 
tor sheet (Fig. 12a). In this case, the mirror image in 
the region z > 0 has no effect on the sheet located at 
z = —d and thus Eq. (43) describes this component in 
the same manner as in the case without a backshort (cf. 
Fig. 8c). 

Slightly more complicated are the second terms in 
Eq. (54), which dominate when p ^$> p. These corre- 
spond to the incident field, E l a , and its mirror conju- 
gate, which corresponds to the reflected field E r a . pass- 
ing through the absorber with negligible interactions 
(Fig. 12b). It is convenient to decompose E r a into two 
parts, one of which, called E r a ,, is phase coherent with 
E l a &tz= —d and satisfies \E r a ,\ = \E r a \ = \E l a \ (Fig. 12c 
and d). Then the entire field is represented with two 
combinations: (c) E\ and E ' a , ; and (d) a superposition 
of K" (= -K’) and E 2 a „ (= E r a ). Note that E\„ \ z= _ d 

has the same phase as E 2 a „ | +d by construction, since 

the image theory forces the field E r a | +d to have a phase 

shifted by n compared to that of E l a |_ = _ . The system 


The contribution from (d) introduces a small violation 
of the symmetry. This component resembles the recep- 
tion pattern of a two-element phased-array antenna with 
uniform (or zero-phase) excitation: 

p | J a" | = cos 0 E | E\, + E\, | = 2 cos 0 E | cos k z d\ \E r \ , 

(56) 

where cos k z d = cos (kd cos 9) is the space factor of the 
phased array. The relation of the current and the inci- 
dent field in the limit of p p (Fig. 12b) corresponds 
to the sum of Eqs. (55) and (56): 

p\J\=2cos9 E \l-e- ik * d cosk z d\\E l \ , (57) 

where the phase factor — arises as a phase differ- 
ence between (c) and (d); the induced current in (c) is 
in-phase with Ei at z = —d while that in (d) is in-phase 
with Ei at z = 0. 

When the surface resistivity is at the neither limits of 
p <C p or p p, the incident electric field corresponds 
to the sum of those in Eqs. (43) and (57): 


2\E t \ 


V 

cos 9h 


P l 

cos Qe 1 — e~ ik * d cos k z d 


|J| , (58) 


which represents the same boundary condition as that 
discussed in Fig. 8 for Eq. (44). Equation (58) is equiva- 
lent to Eq. (54) obtained from the exact calculation done 
in Appendix C. 

When p = 77, the E- and iT-plane responses are maxi- 
mally symmetric and the residual arises from the phased- 
array type contribution of Eq. (56). More strictly, the 
systematics arises due to the difference of this compo- 
nent between E and H planes; the difference of the coef- 
ficient in Eq. (56) is [1 — cos0] cos k z d since 9e = 9 ( 9e = 
0) for .E-plane (iT-plane) incident. This is consistent 
with the fact that S hs (p = 77, d, 9) oc [1 — cos 9] cos k z d as 
in Eq. (47). 

We close the discussion here by pointing out that a 
reflective backshort made of perfect magnetic conductor 
(PMC), as opposed to perfect electric conductor (PEC) 
or metal, would eliminate the polarization systematics. 
With the PMC backshort, the ideal configuration of the 
absorber sheet and the backshort is to place them in- 
finitesimally close to each other. In this case, the phased- 
array type contribution vanishes and the symmetry be- 
tween E and H planes is recovered. Such a configu- 
ration is indeed mentioned in the current sheet model 
as an implementation of an ideal current sheet [48, 49]. 
However, a practical implementation of a PMC tends 
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to have restrictions such as finite bandwidth. This is 
because an implementation of PMC often consists of a 
PEC and a quarter wave delay, and the latter intro- 
duces a wavelength dependence. Indeed, one can regard 
the PEC backshort discussed above as a PMC backshort 
right behind the absorber, where the spatial distance d 
is the quarter wave delay. A phase delay realized by the 
spatial separation depends both on the incident angle 
and on the frequency, and only exactly realizes the de- 
sired boundary at a single frequency. Thus, it leads to 
the systematics discussed above in detail. 

5. Diffraction Effect due to Finite Size of the Ab- 
sorber 

In this section we argue that diffraction effects due to 
the finite size of the absorber are minor. Such diffrac- 
tion effects are discussed in detail in Refs. [50-52], Ref- 
erence [51] shows that the net differential response can 
be approximated as ~ (A/p) • 9 %, where p is the size of 
the absorber. As expected, the larger p is, the smaller 
the differential response. That analysis also shows that 
the differential response has a cos2/> (sin 2<p) pattern 
for a channel sensitive to Q ( U ) polarization. This is 
significant for the purpose of CMB polarization mea- 
surement, as the spurious polarization sourced by this 
effect does not create the so-called R-mode, or parity 
odd, pattern [32, 33]. 

The analysis in Ref. [51] is for a blackbody absorber or 
radiator. However, the difference between an electrically 
thin resistive absorber of interest and a blackbody ab- 
sorber is minor, and only involves angular response pat- 
tern, where the blackbody absorber follows cos 9 while 
the absorber sheet follows a different functional shape, 
Eq. (39), for example. Their difference is small (see, 
e.g., Fig. 6), and is only significant for large angle in- 
cident waves, 9 > 60°, which would either be stopped 
by baffling or have little contribution to the total solid 
angle. Thus, the diffraction estimate for the blackbody 
absorber serves as a good order-of-magnitude estimate 
for an electrically thin resistive absorber, too. A more 
rigorous estimate of the thin absorber response is possi- 
ble with the extension presented in Ref. [52] . 

We also note that these analyses adopt metal bound- 
ary conditions, taking the frame around the absorber as 
a perfect conductor. Strictly, the diffraction effect de- 
pends on the boundary condition and the metal bound- 
ary is one of the worst cases. The diffraction pattern 
in general depends on the impedance of the aperture 
frame [53]. 

There are other possible effects that arise due to the 
finite size of the absorber. For example, the edge of 
the absorber, which may have different impedance from 
the main absorbing area, could affect the image quality. 
Effects of this type are dependent on specific implemen- 
tations of devices and the configurations in which the de- 
vices are placed, such as edge tapering of the absorber re- 
sistivity and baffling around the absorber surface. These 
implementation dependent effects are beyond the scope 
of this paper. 


6. Conclusion 

Estimates of polarization beam patterns are presented 
for a multimoded bolometer employing an electrically 
thin resistive absorber at the limit of infinitely large 
area {N mo & es 1). For a freestanding thin absorber, 
cross polarization and differential response can be elim- 
inated by choosing the surface resistivity to be half of 
the free-space impedance [Eqs. (34) and (35)]. The ab- 
sorption efficiency can be significantly improved by em- 
ploying a reflective backshort termination. Although the 
backshort introduces non-zero cross polarization and dif- 
ferential response, the levels of these systematics can be 
suppressed by choosing the resistivity matched to the 
free-space impedance and optimizing the position of the 
backshort. For a practical application where the absorb- 
ing area is finite, diffraction sources additional beam sys- 
tematics. However, the differential response due to the 
diffraction is small when the number of modes is large. 

The low levels of cross polarization lead to high effi- 
ciency in polarization detection. The smallness of the 
differential response is critical for an accurate measure- 
ment of polarization, as it can create spurious polariza- 
tion from fluctuations in intensity. One significant po- 
tential uses for these sensors is to measure the signature 
of inflation during the early epic of the universe, which 
manifests itself as odd-parity patterns in the CMB po- 
larization [30, 31]. The differential responses due to the 
two types of residual systematics mentioned above, one 
due to a backshort and another due to diffraction, have 
angular patterns of even parity and do not contaminate 
the inflation signature [32, 33]. 
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Appendix A: Explicit Definitions of Unit Vectors 

Unit vectors that define Cartesian coordinates are: 




(Al) 


Unit vectors for polar coordinates are: 


sin 9 cos (f> 
sin 9 sin < 
cos 9 


ee = | 

cos 9 cos (j) 
cos 9 sin 4> 
— sinf? 

— sin (j) 

COS (j) 

1 


(A2) 
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Here are unit vectors that describe the fields of the two The reflected and transmitted fields can similarly be 
polarization bases in Ludwig’s third definition [39]: written as 


e v = eg cos <f> — e<£ sin <p 

( cos 9 cos 2 (j) + sin 2 </> 

— (1 — cos 9) sin cj) cos (j) 

— sin 9 cos </> 

e.h = eg sin <f> + cos cj) 

( —(1 — cos 9) sin ^ cos (j) 
cos 9 sin 2 (j) + cos 2 (f> 

— sin 9 sin cj) 

They are perpendicular to e r by construction and satisfy 
e r x e v = eg, and e r x 4 = — e v . (A5) 




Parallel and perpendicular polarization vectors for x-z 
plane incident waves are: 


(A6) 




They are related to (eg, e $) and (e„, e^) by 

R z (4>)e\\ = eg = cos (f>e v + sin fieh , 
Rz((/>)e± =&</, = — sin (j>e v +cosc/)e h , 


(A7) 


where R z ((f>) is a rotation matrix around the z axis: 
cos (f) — sin (j) 0 

R z (<j)) = | sin^ cos 0 | . (A8) 

0 0 1 


Appendix B: Analytic calculation of the response of 
freestanding thin absorber 

We consider the arrangement shown in Fig. 4. The elec- 
tric and magnetic fields of an arbitrary polarized plane 
wave can be expressed as a linear combination of verti- 
cally and holizontally polarized components with ampli- 
tudes a v = a(9, cj), v ) and ah = a(9, <f>, h ): 

Ei(r; k) = e ~ lk r ( a v e v + a h e h ) 

= e~ lk r R z (cj)) (a.|| e|| + aj_ej_) , 

B i (r-,k) = ie~ ikr (a v e h -a h e v ) (B1) 

c 

= ^e~ lk r Rz(4>) (a||e_L - aj_e||) . 

Amplitudes an and a± are defined for convenience and 
related to a v and ah as 



with 


(B2) 


E r (r- k) = e~ ik r R z (<t> + n) (fo||e M + b ± e ± ) , 
B r (r;k) = -^-e~‘ kr R z ((j) + tt) (fo|| e_i_ -6j_e||) 
and 

E t (r; k ) = e~ lk r R z (<j>) (c|| e|| + ciei) , 
B t (r;k) = l e ~ ik ' r R z (4>) {c\\e±_ - cj_e||) , 


(B4) 


(B5) 


respectively, where k = (k x , k y , — k z ) T ] and fo|| _l and 
C|| i± are the amplitudes of the transmitted and reflected 
waves, respectively. 

The absorber can either be a membrane or a pair of 
orthogonal resistive wire grids running in the x and y 
directions. For the latter, we allow each grid to have 
different resistivity, including the case where one of the 
grids is absent (i.e. , infinite resistivity). At z = 0, where 
the absorber is, all three plane wave components have 
the phase e~ lt ' kxX+kyV ' > and thus the surface current den- 
sity on the absorber induced by the incident field can be 
written as 


J (x, y , z) = S(z)e < kxX + ky v) (j x £ x + j y e y ) , (B6) 

where j x and j y are complex current amplitudes. 

The vector potential A of the field induced by the cur- 
rent J is the solution of the following Helmholtz equa- 
tion: 


V 2 A + k 2 A = -pj , (B7) 

with magnetic permeability y. Note that we implic- 
itly choose the Lorenz gauge by adopting the Helmholtz 
equation. The solution is 

A = (j x e x +j y e y ) , (B8) 

and the corresponding electric and magnetic fields (z ^ 
0) are 

E(r\ k) = ^~y x (V x A) , 

ik (B9) 

B(r; k) = V x A . 

Thus, the field discontinuity at z = 0 due to the current 
J is 


E 


- E 


z- S' — 0 


2 -> + 0 


= ye lkxX lkyV sin 9 ( j x cos (j) + j y sin <j>) e z , 


B 


B 


>-o 


2 -> + 0 


—ik x x—ikvV f • ~ i • \ 

— ( Jy^x Jx^-y) 


(BIO) 


R = 


( COS</> 

l sin (f> 


— sin (j) 

COS 4> 



(B3) 


where y = \fyje is the wave impedance the surrounding 
medium (e.g., vacuum). The electric and magnetic fields 
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of the left hand side are related to Eqs. (Bl), (B4), and 
(B5) evaluated at z = 0. This leads to the following 
boundary conditions on the field amplitudes: 


a,| -6|| -C|| =0, 
a± — b± — c± = 0 , 


(BH) 


and 


“II + h \\ ~ C \\ = V C7a C0S ^ + 3y sin />) , 

Tf 

a± +b± - c± = (j x sin 0 - j y cos 0) , 

cos 9 


(B12) 


where the former [latter] comes from the electric [mag- 
netic] field component of Eq. (BIO). 

In addition to the Maxwell equations used above, the 
current density and the electric fields at * = 0 are related 
by Ohm’s law, leading to 


( Pxjx \ _ rj( C|| cos 6» \ 

V PvJy J V C -L / 


(B13) 


where p x and p v are the surface resistivities of the ab- 
sorber in the x and y directions, respectively. 

Solving Eqs. (Bll), (B12) and (B13), we obtain 


c ll 

c± 



(B14) 


h 

b± 



(B15) 


jx 

jy 



(B16) 


where T , 7Z, and A are transmission, reflection, and 
absorption matrices defined as 


T = 


1 + ^-QR 1 GRQ 


(B17) 


7Z =— sec 6QR GRQT 


2 cos 


9(QR~ 1 GRQ ) 1 


-1 


(B18) 


A. = r]- 1 


Q 1 R~ 1 G~ 1 



with 


Q = 


cos 9 0 \ 

0 1 ) ’ 


G = 


V/Px 0 \ 

o p/py ) ' 


Combining Eqs. (B2) and (B16), we obtain 


(B19) 


(B20) 


jx 

jy 


= AR l 




(B21) 





^0 


Fig. 13. A transmission line of impedance Zq with a re- 
sistance R, modeling a freestanding absorber with an on- 
axis incident wave. The incident wave with a voltage am- 
plitude v induces current j through the resistance, where 
j = v-2/(2R + Z 0 ). 


where 


s -i = 1 ( a+(Px) P(p y ) \ 

4 cos (9 V P(px) a-(py) J ’ 


(B22) 


with 


a±(p) = (2 p + rj)( 1 + cos 6) ± cos2^)(2p — rf)( 1 — cosf?) , 
(3(p) = sin 20(2,0 — 77) (1 — cos 9) . 

(B23) 


When the matrix S is diagonal, the vertical (horizontal) 
polarization amplitude a v (ah) only couples to the x ( y ) 
current amplitude j x ( j y ) and thus cross polarization 
is zero. Further, when the two diagonal elements of S 
are equal, the angular response patterns to vertical and 
horizontal polarization waves are the same and thus the 
differential response is zero. 

For on-axis incidence with 9 = 0, the relation between 
the incident field and the current is equivalent to the 
case for a one-dimensional transmission line (Fig. 13), 
as expected: 

X — ( ^/^P x 0 \ ('B241 

S *=o-V 0 2/(2 p y + v)J- (B24) 


Appendix C: Analytic calculation of the response of 
thin absorber with reflective backshort termination 

We adopt image theory and use the setup of Fig. 9b. As 
the mirror conjugate flips the sign of the z components 
of vectors and the sign of the electric field, the conjugate 
field of Ei can be written as 


Ej = -e lk r (a v e.v + a h e.yj) 

= e~ tkr R z (</) + ir)(a\\e\\+a±e ± ) , 

with 


(Cl) 



&h = 



e h . (C2) 


The conjugates of the reflected and transmitted compo- 
nents, Ef and Ef, respectively, and all the associated 
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magnetic fields can be defined in the same way. In par- 
ticular, those for the transmitted field are 

Ef{r ; k) = e~ lkr R z {<j> + tt) (c||c§|| + c±e±) , 

1 ■- (C3) 

Bf(r;k) = --e lkr R z (<f> + 7r) {c\\e±_ - c_Ley) . 

The boundary condition is defined in the same manner 
as Eq. (BIO) but at the surface of z = —d, leading to 


7 d|| - 7 * 6 || - 7 C|| + 7 *cy = 0 , 
7aj_ — 7*6 _l — 7c_l + 7*c_l = 0 , 


(C4) 


and 


7 a || + 7*^|| - 7 C || - 7 *C|| = r)(j x cos 0 + j y sin 0 ) , 

fj 

7 a _L + 7*fr-L - 7 C -L - 7 *c_l = (jx sin 0 - j y cos 0) , 

COS u 


(C5) 


with 7 = e Ite * . On the other hand, Ohm’s law leads to 


Pxjx 

Pyjy 


C || cost 
C_L 


= (7 — 7 *)R 
Solving Eqs. (C4), (C5) and (C 6 ), we obtain 
= "Abs 


a \\ 

a± 


(C 6 ) 


(C7) 


with 

A. bs = P _1 


1 




7 — 7 ’ 27 

Combining this with Eq. (B2), we obtain 




= «4.bs-R 1 


a„ 

a/i 


= £ 


bs 


a„ 

a/i 


where 


~—ik z d 


= , C^( Px ) p hS (Py) 

4sinfc z <icos0 \P hs (Px) a bs (Py) 


(C 8 ) 

(C9) 

(CIO) 


£ _1 
with 

a^ s (p) = [(p + p) sin k z d — ip cos k z d ] (1 + cos 0 ) 

±5 hs (p,d,6) cos 20 , (Cll) 

/3 bs (p) =S hs (p,d,6) sin 20 , 
and 

<5 bs (p, d, 9) = [(p — 7 ) sinfc^d — *pcosfc z d] (1 — cosf?) . 

(C12) 

Again, for on-axis incidence, 0 = 0, the relation be- 
tween the incident field and the current is equivalent to 
that of the one-dimensional transmission line (Fig. 14), 
as expected: 


J bs 


with 


fo (p) = 


co(Px) 0 
?=0 V 0 C 0 (p y ) 


2e ikd sin kd 

(p + rj) sin kd — ip cos kd 


(CIS) 


(C14) 




d-y 4 


Fig. 14. A transmission line of impedance Zq with a re- 
sistance R and a reflective termination, modeling an ab- 
sorber with reflective backshort termination and with an 
on-axis incident wave. The distance between the resis- 
tance and termination is d and the wavelength in the trans- 
mission line is A. The incident wave with a voltage am- 
plitude v induces current j through the resistance, where 
j = v ■ 2e zkd sin kd/[(R + Zq) sin kd — iRcoskd], with wave 
number k = 2n/X. 


Appendix D: Possible Systematics due to the Finite 
Physical Size of the Grids 

All the discussions in this paper are based on an assump- 
tion that the absorber can be modeled as a thin resistive 
membrane. When the absorber is a pair of orthogonal 
grids, it deviates from the ideal membrane model in the 
following three ways: 1 ) there is a finite distance between 
the layers of orthogonal grids, 2 ) the grids have non-zero 
pitch and consist of wires with non-zero cross-sections, 
and 3) the two orthogonal grids may couple to each other 
through near-field effects. In this Appendix, we estimate 
the magnitude of the systematics due to these aspects 
and discuss the conditions on the physical dimensions of 
the device required to suppress such artifacts. 

1. Non-zero Distance between the Two Layers 

When the absorber is a pair of resistive grids, which 
are sensitive to orthogonal polarizations, there is a fi- 
nite gap between the two grid layers. This may lead to 
additional beam systematics of cross polarization and 
differential response. To estimate the effect, we consider 
a freestanding absorber with p = p/ 2 , where there is no 
systematics in the nominal configuration. 

We first point out that the primary effect here is the 
cross polarization, not the differential response. This 
can be seen in Eq. (32); when p x = p/2, the response of 
the grid running in the x direction, G v COl is independent 
of p v including the case where the grid running in the y 
direction is absent, or p y = 00 . Thus, the co-polar beam 
shape does not see the effect of the other grid to first 
order. 

Cross-polarization, on the other hand, can arise from 
the finite gap. Consider the electric fields of incident, 
reflected, and transmitted waves projected on the ab- 
sorber surface plane (x-y plane), E t | , E r | , and E t | 
respectively, for vertical polarization incident. Fig- 
ure 15a shows an example with an incident angle of 
0 = 45°, where the cross polarization is maximum, 
and 9 = 50°. The boundary condition guarantees 
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E i | + E r | = E t | . Zero cross-polarization of the nom- 
inal configuration corresponds to a vanishing y compo- 
nent of the field: Ei tV + E r , y = E ttV = 0. However, this 
cancellation of the y components of the incident and 
reflected fields is only exact on the surface of the grid 
running in the x direction, which we define as z = 0. 
When the other grid running in y direction is placed at 
z = —e (Fig. 15b), the y grid feels a residual electric 
field of order 

£ 7T6 

(1 — cos 6) sin k z e ~ 27 r— cos 9(1 — cos 9) < — . (Dl) 
A ZX 

Thus, the cross-polar level is 

XP ~ (7re/2A) 2 . (D2) 




Fig. 15. (a) The electric field of a vertically ( x ) polarized in- 

cident wave with ( 9 , <f>) = (50°, 45°) projected on the plane of 
the absorber (the x-y plane), (b) A schematic figure showing 
the case with a gap e in the z direction between the planes 
of x and y grids, with the vertically polarized incident wave. 


2. Finite Pitch and Cross Section of the Grid 

Compared to an ideal resistive sheet, a grid of resistive 
wires has non-zero spacing between the grid wires and 
non-zero cross section of each wire. These can cause de- 
viations from the ideal current sheet when the physical 
dimensions are not small compared to the wavelength 
A. The requirements for the smallness are discussed in 
detail in literature (see, e.g., Refs. [54, 55] and refer- 
ences therein) in the context of grids made of conductive 
wires. Estimates for the grid of resistive strips do not 
significantly deviate from conductive wires when the re- 
sistivity is similar to or less than the impedance of free 
space, p < 77 [56, 57], which is the parameter space of 
interest here. In the regime X^$> g > Ana, where g is the 
spacing between the wires and a is the radius of each 
wire of circular cross-section, the first order estimates 
of the cross-polar level and the differential response are 
(see, e.g., Ref. [54]) 


/ 27 r 2 a 2 A 

\~W) ’ 


DR 
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A 2na ) 
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An estimate for a grid made of thin strips of width w 
can be obtained by substituting w = 2a, to first or- 
der. In the regime A 3> 47ra > g , the approximation of 
the term involving lngf/27ra becomes invalid as pointed 
out in Ref. [55]. However, the numerical analysis in the 
reference shows that the term is still a monotonically de- 
creasing function of a/g. Thus, one can still put a rough 
upper limit on the systematics based on Eq. (D3) and the 
monotonic dependence on a/g even in this regime. In the 
implementation of Ref. [26], for example, the parameters 
are a ~ 1 /iin and g ~ 30 gm, and satisfy A > j > 47ra. 
The systematics in this case are XP ~ (l/\[gm\) 2 and 
DR ~ (10/A[/xm]) 2 , negligibly small for millimeter and 
submillimeter wavelengths. 


3. Near-field Coupling 

All the discussions above assume only the far-held ef- 
fects, or the radiating modes, and we ignored the near- 
held effects due to the evanescent modes. Here we show 
that the near-held coupling between the orthogonal grids 
is small when the distance between the grids, £, is sim- 
ilar to or larger than the grid spacing, g. Owing to the 
periodic symmetry of the system, the scattered held in- 
cluding the evanescent modes can be expanded in terms 
of a Floquet series [57, 58]. The phase component of the 
Floquet series is 

g-jfc z ,mn |Az| g — X+k y , n y) (D4) 


with 


k x ,m = k x + 2-Km/g , k V}n = k y + 2i rn/g , 

kz,mn = y k 2 k x ,m k y ^ n (Im[/c z ,mn] £ 0) 5 


(D5) 


where A z is the distance from the scattering grid, and 
m and n are the integer numbers corresponding to the 
series indices. Our interest here is the non-radiating 
modes, or the modes with m / 0 or n / 0. When 
the grid spacing is small compared to the wavelength 
(g • C A), the non-radiating mode has 


k z ,mn - 
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(D6) 


and the held strength decays as e~ 2 ^^ Az ^ 9 or faster. 
Thus, the near-held coupling between the paired grids 
is suppressed by a factor of e -2 ’ 1 ' 6 / 9 . This factor is small 
when £ is similar to or larger than the grid pitch (e > g) . 
Therefore, the near-held effect can safely be ignored in 
this regime of e > g, with g < CA. 

Appendix E: Relation among Resistive Sheet, Conic 
Reflector, and Parallel Current on Their Surface 

Koffman showed the condition where the induced cur- 
rent on a conic rehector illuminated by a feed hows in 
parallel [47]. The condition is met when the magnetic 
held H at the rehector surface satishes 

Hn re + cos 6 

= tan cp ■ 


XP 


1 + re cos 9 


(El) 
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where the angle (6, cf>) specifies a position on the reflec- 
tor with the focus as the origin of the coordinates, H\\ 
and H± are the magnetic field along eg and ey,, respec- 
tively, and k is the eccentricity specifying the geome- 
try of the conic reflector: k = 0, 0 < k < 1, n = 1, 
1 < k < oo, and k = oo correspond to sphere, ellipsoid, 
paraboloid, hyperboloid, and plane, respectively. Koff- 
man pointed out that a feed radiation pattern satisfies 
Eq. (El) when it can be expressed as a superposition 
of electric- and magnetic-dipole radiations, and relative 
strengths between them coincide with the eccentricity n. 
Note that the object is at infinitely far for astronomical 
applications and thus a paraboloid reflector or an equiv- 
alent system [59, 60] is usually employed. This is why 
we adopted Ludwig’s third definition, in which a pure 
polarization pattern satisfies Eq. (El) with k = 1. 

We can see two connections between this symmetry 
explored by Koffman and our result for a freestand- 
ing resistive absorber sheet. One of them can be seen 
by replacing the feed in Koffman’s setup by a resistive 
sheet. In the replacement, we also replace the radia- 
tion field pattern from a pure-polarization feed by the 
absorption field pattern that induces current in only 
x, i.e., j y = 0. Equations (A7), (Bl), and (B2) lead 
to --ff||) T = (a||,a_i_) T = R(a v ,a h ) T . Thus, 

Eqs. (B21), (B22), (B23), and j y = 0 lead to 


Hu 


H ± 


tan <f> ■ 


J- P + C0S 6 
1 + #■ cos 9 


(E2) 


for the reception pattern of the thin resistive absorber. 
Equations (El) and (E2) have equivalent form, due to 
the following parallel between Koffman’s and our results: 
k — > 0 (p — > oo) corresponds to magnetic dipole (mag- 
netic conductor sheet), while k — > oo (p — > 0) corre- 
sponds to electric dipole (electric conductor sheet); the 
electric and magnetic contributions balances when k = 1 
( p = rj/ 2). See also the discussion at the end of Sec. 4. A. 

Another connection between Koffman’s and our re- 
sults can be seen by replacing the conic reflector in 
Koffman’s setup with a freestanding thin resistive ab- 
sorber. Namely, one illuminates the resistive absorber 
by a superimposed electric- and magnetic-dipole radi- 
ations with their relative strength of k. According to 
the derivation above, the induced current on the ab- 
sorber flows in parallel when the sheet resistivity sat- 
isfies k = r]/2p. Thus, instead of the eccentricity of 
the conic reflector, one can use the resistivity of the ab- 
sorbing sheet as the parameter tuned to match with an 
arbitrary superposition of the dipoles. 
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